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I.  INTRODUCTION 


Measured  dipole  element  patterns  in  a  number  of  phased  arrays 
indicate  the  appearance  of  blindspots  [1,2].  Analyses  which  disre¬ 
gard  the  presence  of  dipole  supports  [3,4]  do  not  predict  such  re¬ 
sonances.  In  order  to  understand  the  cause  of  these  blindspots  and 
how  to  avoid  them  it  is  necessary  to  determine  the  dipole  element 
performance  with  inclusion  of  the  feed.  To  circumvent  the  ensuing 
complication  of  simultaneously  satisfying  boundary  conditions  on  the 
dipole  arms  and  on  their  supports,  it  is  desirable  to  first  determine 
the  modes  of  the  feed  region. +  These  modes  individually  satisfy  the 
boundary  conditions  on  the  supports,  as  well  as  the  Floquet  require¬ 
ments  on  the  unit  cell  walls.  The  analysis  of  a  dipole  array  scan 
performance  will  then  be  reduced  to  an  evaluation  of  the  discontinuity 
presented  by  dipole  arms  between  two  different  unit  cell  waveguides 
(one  representing  the  feed  region  and  the  other  the  air  region  above 
the  array). 

Section  II  is  devoted  to  the  derivation  of  the  dispersion  relation 
and  the  fields  for  the  TM  modes  of  the  feed  region  of  a  dipole 
phased  array,  in  strip-line  geometry  and  with  a  rectangular  lattice. 
Section  III  discusses  the  TEM  modes  supported  in  this  feed  region. 
Section  IV  addresses  the  feed  region  modes  for  the  case  of  a  triangu¬ 
lar  array  grid.  Section  V  discusses  details  of  the  numerical  analysis. 
Section  VI  considers  the  numerical  results  and  presents  a  number  of 


+  A  similar  approach  was  taken  in  ft.  Lewis,  A.  Hessel,  G.  Kmttel, 
"Performance  of  a  protruding-dielectric  waveguide  element  in  a  phased 
array,"  IEEE  Trans.  Antennas  Propagat.,  Vol.  AP-20,  pp.  712-722, 
Nov.  1972. 


-1- 


conclusions.  The  appendix  is  devoted  to  the  orthogonality  of  the 
degenerate,  feed  region  TEM  modes. 


I' 

I 

I 


II.  TMZ  modes  of  the  feed  region  for  a  rectangular 

LATTICE 

A  typical  dipole  array  element  is  shown  in  Fig.  1,  and  a  simpli¬ 
fied  ribbon  model  of  the  element  is  seen  in  Figure  2. 

The  feed  region  part  of  the  unit  cell  may  be  regarded  as  a  finite 
section  of  the  infinitely  extended,  uniform  guiding  structure  pictured 
schematically  in  Figure  3. 

The  strip  width  is  assumed  small  compared  to  the  free  space 
wavelength.  Therefore,  the  usual  approximation  of  only  longitudinal 
current  is  made.  This  in  turn  implies  that  only  the  TMZ  unit  cell 
modes  are  affected  by  the  strip  loading.  Hence,  only  the  TMz  propa¬ 
gation  characteristics  and  the  associated  modal  field  distributions  need 
be  determined.  The  strips  are  invisible  to  the  TEZ  fields  which  are, 
therefore,  represented  in  terms  of  the  usual  Floquet  unit  cell  modes. 

To  determine  the  TM  modes  two  coupled,  linear,  homogeneous 
integral  equations  are  set  up.  These  equations  express  the  Floquet 
boundary  conditions  on  the  unit-cell  walls  as  well  as  the  vanishing  of 
the  tangential  E  field  on  the  conducting  strips.  They  are  solved  by 
Galerkin's  procedure,  resulting  in  modal  fields  that  exactly  satisfy  the 
Floquet  condition  but  only  approximately  fulfil  the  requirements  on 
the  strips. 

To  set  up  the  integral  equations,  the  axial  component  of  the 
electric  field  is  expressed  in  terms  of  the  Green's  Function  for  an 
infinite,  electric-current,  line  source  located  in  the  unit  cell-wave- 
guide,  with  e^wt‘KZ^  variation  and  the  allowed  values  of  k  being  as 
yet  undetermined .  Accordingly , 


(1) 


E  (r;  k,  kto)  =  e-jKZ  [f  G(£  £*;  k.  ktQ)  JjCfi’)  dp'  + 
z  "  LU  strip  1  — 

S  G(e.e’;  K.  kt0)  J2  (£')  dp'  ]  . 
strip  2  — 

In  (1)  the  strip  current  densities  are  2(2') »  an(*  un**  ce^ 
Green's  function  is 


=  ^  JAn  Y"  (X'X')e 


,JVy-y,) 


2nn 


where  ^  =  V  +  ^ 

kxn  =  kkn  5  0 

"Jkxn(x'x,)  „  jkxn(x-x,)  -ik  lx-x'| 


(2) 

(3) 

(4) 


Y»(x'x'>  -  iosc-^7  *  7*^7;  * e 

c  “i  e 

,  W 

n  =  kyo  =  k  sin  6Q  sin  *Q 

(6) 

£  =  k  sin  0Q  cos  4>0 

(7) 

Ko  =  A  +  A 

(8) 

£  =  x  xQ  +  y  y0  ,  £'  =  x'  Xq  +  y' 

(9) 

r  =  £  +  z  zQ 

(10) 

e  is  the  free  space  permitivity . 
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The  angles  (0Q,  $Q)  denote  the  usual  beam  pointing  angles  of  a 
phased  array  whose  broadside  direction  is  zQ.  It  should  be  stressed 
that  k  d  is  not  the  steering  phase  shift  along  x,  which  is  given  by 
£d  •  rather,  k  and  the  other  k  are  propagation  constants  for 

A  Aw  All 

E-type  modes  of  a  phase-shift-wall  parallel-plate  waveguide  of  height 
dy.  These  modes  travel  in  the  x-direction  and  have 

-j(k^ny+<2) 

e  1  transverse  variation. 

By  enforcing  the  boundary  conditions  of  the  strips,  i.e. 


Ez(+  Sf.O)  =  0  for  |y|  <  \ 


the  following  coupled  integral  equations  for  2  are  obtained  from 
(1): 


G(-J,  y,  -f .  y';  •c.k^Cy')  +  G(-£,  y,£,  y‘;  K,ktQ)  J2(y')  =  0 


for  |y|  <  % 


Gf  ,  y,  -f.y';  k,  k^)  J^y')  +  G(£.  y,^  ,  y’;  k,  k^)  J2(y')  =  0 

h 

where  GJ  a  J  G  J  dy'  (see  Fig.  3).  (13) 

~  h 

"  I 


In  general,  to  solve  (12)  and  (13),  the  current  densities  would 
be  expanded  in  a  set  of  linearly  independent  basis  functions.  How¬ 
ever,  in  view  of  the  narrow  strip  width  only  one  term  of  the  expan¬ 
sion  is  used  here,  i.e., 
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Ji(2(y')  =  or1)2<i>(y')- 


(14) 


Application  of  Galerkin's  method  now  yields  the  equations 

^Z11(K'  W  212(k'  ktoA  /“] 
i221(K'  ktQ)  Z22(k,  ktQ)  J  U2 


=  0 


(15) 


where  jZ  (k,  ktQ)  =  <Ky),  G((-1X  ^  ),  y,  (-l? (^  ),  y';  k,  ktQ)  <l»(y')>. 


h 

2 


and  <f,g>  =  J  fg  dy 
h 

■  2 


(16) 

(17) 


(the  form  jZpq  is  introduced  so  that  Z^  and  Z22  are  real). 


In  view  of  the  following  symmetries  in  (5): 


Yn(x\x')  =  yn(-x',-x') 


(18) 


and 


Yn(x,x') 


'rtn(x',x)' 


(19) 


equations  (15)  can  be  written 


'zn(K,  ktQ)  212(k,  ktQ)' 


^i2(k,  ktQ)  Z11(k,  k tQ)y  \a2 


(20) 


For  the  choice  ijj(y')  =  1,  i.e.,  a  constant  current  distribution  on 
each  strip, 


»  1 


sink*A 


2u(k,  V  ^  •  coskxndx-cosi;d 


ft sin  kyn11  \ 

5? 


/ 


(21) 
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2 


»  _1_  sink  (d-w)+ej^x  sink  w 
Z12(k'  V  kxn  c°sl£x_ndx-c0s4dx 


'sinkynh 

■ 


Alternatively,  if  4i(y ')  is  chosen  to  be  the  static  distribution 


>Ky')  = 


/c£)2  -  y2 


which  satisfies  the  edge  condition  on  the  strip,  then  the  zeroth  order 

/  kynh v 

Bessel  function  of  first  kind  JQ  (¥)  ,  replaces  — J  in 

and  (22). 

Scan  dependent  solutions  for  2  exist  only  for  values  of 
k  satisfying 

211(K'  V  "  |Z12(K'  V|2  =  0  '  (24) 

which  is  the  desired  dispersion  relation.  Its  solutions  yield,  to 
within  our  approximations,  the  propagation  constants  of  the  TM 
modes,  guided  by  the  two-conductor  strip-line  in  a  unit  cell.  The  k’s 
depend  on  frequency,  on  geometry  and,  unlike  in  closed  waveguides, 
on  scan  variables.  It  can  be  shown  in  the  usual  manner  that  the 
values  of  k  are  either  real  or  purely  imaginary. 

For  a  particular  k,  the  associated  2  are  evaluated  from  (20) 
to  within  a  normalization  constant,  and  the  relative  current  distribu¬ 
tion  on  the  strips  is  thus  determined.  Subsequently,  E  is  calculated 
from  (1)  and  other  field  components  are  obtained  from  Maxwell’s 
equations.  The  modes  so  derived  are  mutually  orthogonal,  although 
they  only  approximately  satisfy  the  boundary  conditions  on  the  strips 
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[5].  It  can  also  be  shown  that  these  TM  modes  are  orthogonal  to 
the  TE2  Floquet  modes  as  well  as  to  the  TEM  modes  derived  in  the 
next  section. 
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III.  TEM  MODES  OF  THE  FEED  REGION  FOR  A  RECTANGULAR 
LATTICE 

The  TEM  modal  fields  cannot  be  determined  in  the  same  manner 
as  the  TM  fields  because  now  Ez  is  identically  zero.  Instead,  the 
electrostatic  potentials  on  the  metallic  strips  are  assigned  constant 
values  ^  and  respectively.  Two  coupled,  inhomogeneous,  linear, 
integral  equations  of  the  first  kind  are  set  up  which  express  this 
requirement.  Application  of  Galerkin's  procedure  yields  a  set  of 
inhomogenous ,  linear,  algebraic  equations  which  are  solved  for  the 
charge  density  n,  or  the  TEM  current  density  J,  where  J  =  ^  q-  To 
formulate  the  integral  equations,  the  static  Green's  function  for  an 
infinite  z-directed  line  charge  in  a  unit  cell  is  employed.  The  Green's 
function  is  given  by 

G(fi.  B':  V  *  5^  ni.^nYn(x’x’)e  ',kyn(y‘y,>  (25) 

with 

-kyn(x-x')  kyn(x-x') 

vx'x’>  -  •  <26> 


The  wavenumbers  kyn,  4  and  ktQ  are  defined  in  (6-8).  If  the  strip 
charge  densities  are  designated  n^(y')  and  q2 (y')»  respectively,  the 
resulting  potential  in  the  unit  cell  is 


*(x,y)  *  J  G(e,  ktQ)n1(y’)dy' 
strip  1  —  1 


♦  J 

strip  2 


G<e.  e';  kt0)n2(y')dy' 


(27) 


r9 


The  boundary  conditions  on  the  conductors  yield  the  integral  equa¬ 
tions, 


♦i  =  g(-  y,  -  y';  kt0)n1(y')  +  G(-^,  y,  y';  kt0)n2(y') 


for  |y  |  <  | 


^2  =  ¥•  ”*2'  y'  Vni(y’)  +  y,  y'; 

where  Gr|  is  defined  as  in  (13). 

Again,  assuming  a  single  mode  approximation 

and  using  the  symmetries  of  (26), 

Y^x'.x')  =  Yn(-x',x’) 


Yn(x,x')  =  Yn(x',x) 
the  Galerkin  procedure  yields 

(P11<V  P12(kto)\  /Ql\  /*l\ 

^•P12(kto)  Pll<kto>/  \Q2^  ^2 


where  <«|»(y),  1>2  Plq(k^)  =  <«Ky).  G(-j,  y,  (-l)q(^),  y’;  kt0)«Ky')> 

(33) 


Ql,2  =  al,2<*<y)'1> 

are  the  electric  charges  per  unit  length  z  on  the  strips. 


For  <Ky')  =  1, 


=  -  1-  I 


sinh  k^,,  dx 


11  -  -  £0  n=.„  2dykyn  cosidx-coShkync 
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L  »  __i_  elCdxsinhkvnw>sinhRvn(dx-w)  /Slnffi 

12t-SS>  ‘o  n—  2dykyn  cos^-coshk^d'  \~Epr 

\  2 

(36) 

In  general,  matrix  E  in  (32)  is  non-singular  and,  therefore,  one  has 
from  (32) 


(37) 


where  the  capacitance  matrix,  C,  is  E  • 


0i 

Choosing  a  particular  combination  (  )  defines  a  mode  since  the 

2 

modal  fields  are  then  determined  via  (37),  (34),  (29),  and  the 
gradient  of  (27). 

In  general,  there  are  two  linearly  independent  choices  for  the 

vector  (^  ),  and  these  yield  two  linearly  independent,  charge  distri- 
2 

butions  for  the  same  TEM  propagation  constant,  k  =  k,  i.e.  there 
are  two  degenerate  TEM  modes.  (This  result  also  holds  when 
more  than  one  term  is  used  in  the  strip-current  expansion). 

A  convenient  basis  for  the  vector  space  {(  )},  is  one  whose 

h 

associated  fields  are  orthogonal  over  the  unit  cell  cross-section.  It  is 
shown  in  Appendix  I  that  the  eigenvectors  of  the  capacitance  matrix 
exhibit  this  property.  The  eigenvectors  are 


^a 


(38) 


-11- 


T 

They  satisfy  <£a  *  $5  =  6a  £,  where  5a  ^  is  the  Kronecker  delta  func¬ 
tion.  The  associated  (real)  eigenvalues  are 

Yb ■  CP11 ;  |pi2l>/4  •  «9> 

where 

A  =  pii  -  |P12I2.  (40) 

As  seen  from  (38),  for  each  eigenvector  the  component  ratio,  which 
represents  the  relative  strip  charge  distribution  for  a  given  TEM 
mode,  is  of  unit  magnitude.  Furthermore,  the  two  modes  have  rela¬ 
tive  charge  distributions  which  differ  in  phase  by  180°. 
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IV.  TM,  MODES  OF  THE  FEED  REGION  FOR  A  TRIANGULAR 
LATTICE 

The  analysis  of  Section  II  is  now  extended  to  the  case  of  a 
triangular  grid  by  choosing  the  unit  cell  shown  in  Figure  4.  The 
Green's  function  (2)  is  also  valid  in  this  case.  By  imposing  boundary 
conditions  on  all  the  strips  a  set  of  four  integral  equations  is  ob¬ 
tained,  which  are  reduced  to  a  set  of  linear,  homogeneous  algebraic 
equations  via  the  Galerkin  procedure.  Imposing  the  steering  phase 
shift  between  the  two  strip  transmission  lines,  shown  in  Figure  4, 
reduces  the  problem  to  the  solution  of  two  equations  in  two  un¬ 
knowns,  which  in  turn  yields  the  dispersion  relation  and  the  modal 
current  distribution. 

In  detail,  the  modal  axial  electric  field  generated  by  the  strip 

A 

current  distribution,  Jj(y')  zQ  (i  =  1,2, 3, 4),  and  with  G  given  in 
(2),  is 

-iKZ  4 

E-(r;K,kto)  =  e  JKZ  1  [J  G(e,e';K,kto)L(y')dy']  (41) 
—  i=l  strip  i  — 

The  integral  equations  are  obtained  by  requiring  E  to  vanish  on 
each  of  the  four  conductors.  In  reducing  these  to  algebraic  equa¬ 
tions  one  term  is  used  in  the  current  expansion,  as  previously  in 
(14). 

For  the  choice  'Ky')  =  1#  the  set  of  equations  is 

(42) 
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i* 


where  the  independent  coefficients  of  the  Hermitian  matrix  are: 


Zll<K'kto>  = 


aa 

I 

n=-«* 


z12^K/kto^  ~ 


oo 

I 

n=-<» 


sin(kxndx)  /sin  ^4 

'  coskxndx*cos^dx  V  V 

V  2 

)Sd  /-jn’S'n11  \2 

sinkxn(dx~w)te  sinkxnw  f _ 2 _  j 

“skxndx-cos«dx  /M4> 


00 

Z13(K'kto)  “  1 
—  n — ® 


*dx  *Vin  kxndx 
e  "2~  -cos-g —  sin  — g - 


coskxndx-cos«dx 


Z14(K'kto) 

—  n — °° 


k  h  ^ 

1  sinkxn(dx/2'wJ+^dxsinkxn(dx/2+w)/^irT~  \ejkvndv 

^  coskxndx'cos^dx  ) 

X_2~  (46) 


J23 


(K'kto>  =  *  k” 
—  n=-«  xn 


sinkxj1(dx/2*w)+^dxsjrkxn(dx/2-w)/sii> 
coskxndx-cos^dx 


e^y 

(47) 


Imposition  of  the  steering  phase  requirements  between  the  two 

ndv 

"K£dx  + 

feed  lines,  i.e.,  of.+2  =  a.e  ,  where  £  and  r|  are  defined  in 

(6-7),  yields 


(48) 
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with 


V-V  =  +  213(~K,kto^ 

22(K,kto)  =  212(K,kt0)  +  214(K,kt0) 

Z3^K,kto^  =  2i2^K,kto->  +  Z23<«'V 

Using  (43-47)  and  (49-51)  one  finds 

VK'V  =  zi  ^K,kto^ 

22(K,ktQ)  -  22(K,kto)  , 
and  (48)  reduces  to 


that 


yielding  the  dispersion  relation  for  the  triangular  grid  case 


Z1  <K'kto>  “  l22(K'kt0)|2  =  °- 


(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 


The  TEM  mode  charge  distributions  for  the  triangular  grid  can 
be  determined  by  a  corresponding  extension  of  the  analysis  presented 
in  Section  III. 

Four  inhomogeneous  equations  are  obtained  via  Galerkin's  proce¬ 
dure  and,  by  imposing  the  steering  phase  requirements,  are  reduced 
to  the  form  of  (32),  where 


pll<kto>  =  IcTT  *  r~ 
y  o  n=-» 


sinhkyn 


d  +2cos(nn)coi 


coshkyndx-cos£d 


(56) 
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P12(kto)  “ 


<*  1 

2 HT  1  k~ 

y  o  n=-®  yn 


sin 


hkvn 


(d  -w)  +  e 


j 


x  _ 

coshkyndx-cos^dx 


sinhkynw 


.W  . £dx 

j  ~2 —  dx  -j  ~2 —  d 

cos(nn)  e  sinhkyn(-2+w)+e  sinhkyn(-£  -w) 

cosh  kyndx-cos£dx 


(57) 


The  charge  distributions  and  fields  may  then  be  calculated  as 
indicated  in  Section  II. 
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V.  NUMERICAL  ANALYSIS 


The  TM  propagation  constants,  which  are  the  roots  of  (24),  were 

evaluated  by  using  a  combination  of  the  step  searching  method  and 

the  secant  method.  For  small  values  of  h/\  (Fig.  3),  convergence 

of  the  series  (21)  was  accelerated  using  the  closed  form  expression 
.  2 

[6]  for  ?  sirc  v-  .  Series  (22)  converges  well  because  its  terms 
n=l  n 

decay  exponentially  with  n.  The  curves  of  Figs.  5-8,  were  obtained 
with  n  ranging  from  -20  to  20.  Extending  the  range  of  n  from  -200 
to  200  did  not  change  the  third  significant  digit  of  either  (21),  (22) 


In  computing  the  curves  shown  in  Figs.  5-8  a  single  constant 
term  was  used  for  the  current  expansion  on  each  strip.  Curves  were 
also  computed  using  the  single  static  term  (23)  and  a  comparison  of 
results  for  these  two  current  choices  is  shown  in  Fig.  9. 

For  the  scan  plane  <|>  =  90°,  the  unit  cell  of  Fig.  3  can  be  bi¬ 
sected  (Fig.  10).  A  magnetic  wall  is  used  because  for  the  range  of 
parameters  over  which  £  was  calculated,  the  currents  on  the  two 
strips  are  in  phase.  This  relationship  is  obtained  by  evaluating  (21) 
and  (22)  (both  are  real  for  $  =  90°)  and  solving  for  2  *n  (20). 
The  bisected  cell  is  identical  to  the  unit  cell  of  a  one-dimensional 
strip  grating  bordered  by  magnetic  walls.  A  circuit  representation  of 
the  grating  is  used  to  construct  an  equivalent  network  [6],  whose 
transverse  resonance  solution  yields  a  dispersion  relation  for  our 
bisected  unit  cell.  The  propagation  constants  obtained  from  this  re¬ 
lation  are  in  close  agreement  with  those  calculated  from  (24),  except 
near  cutoff. 
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The  phase  angle  of  the  relative  current  distribution  vs.  scan  is 
shown  in  Fig.  11  for  the  dominant,  plus-minus,  TEM  mode.  The 
phase  angle  for  the  plus-plus  TEM  mode  is  obtained  from  Fig.  11  by 
adding  180°  to  each  phase  value  shown.  The  angles  are  measured 
using  the  left  strip  (Fig.  3)  as  the  reference.  The  relative  charge 
distributions  were  calculated  using  eigenvectors  (38).  Series  (36) 
was  evaluated  with  n  ranging  from  -20  to  20. 

Exactly  at  broadside  the  TEM  plus-plus  mode  does  not  exist, 
since  all  the  conductors  are  at  the  same  potential. 

For  the  triangular  grid,  the  TM  propagation  constants  (Fig.  12) 
and  the  TEM  current  distribution  (Fig.  13)  were  calculated  in  similar 
fashion  except  that  (21)  was  replaced  by  (49)  and  (22)  by  (50).  The 
series  of  (43),  contained  in  (49),  is  identical  to  that  in  (21)  and  its 
convergence  can  be  accelerated  as  discussed  above.  The  other  series 
in  (49)  and  (50)  all  converge  rapidly. 

The  calculations  for  the  rectangular  and  triangular  grids  were 
performed  by  computer  programs  which  were  checked  against  each 
other  for  the  special  case  of  the  strip  separation,  w,  equal  to  one- 
half  the  unit  cell  x-dimension,  since  in  that  case  the  two  grids  are 
identical . 
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VI.  NUMERICAL  RESULTS  AND  CONCLUSIONS 

Figures  (5-8  and  12)  indicate  that  under  appropriate  conditions  a 
propagating  TM  mode  exists  in  the  feed  region  of  a  dipole  array. 
The  normalized  propagation  constant  of  this  mode  for  a  rectangular 
grid  geometry  is  shown  in  Fig.  5-8  for  different  values  of  d  A, 
d  A,  wA  and  hA,  respectively.  It  is  seen  that  £  increases  when 
the  unit  cell  dimensions  increase,  or  when  the  strip  separation  wA, 
or  the  strip  width  hA,  decrease.  Figures  5-8  show  that  £  decreases 
monotonically  with  azimuthal  angle  0  and  elevation  angle  6.  However, 
in  Fig.  6,  £  decreases  monotonically  with  elevation  angle,  only. 

In  Fig.  12,  £  is  shown  vs.  scan  for  two  equilateral  triangular 
grids.  Here  too,  £  increases  as  the  unit  cell  dimensions  increase. 

The  TM  strip  currents  are  not  shown  in  the  figures.  However, 
knowledge  of  the  modal  propagation  constants  makes  it  possible  to 
evaluate  the  currents  from  (20).  For  each  mode,  (24)  must  be  satis¬ 
fied  and  2^  =  ±  |Z12I-  By  substituting  this  result  into  (20)  it  is 
seen  that  the  modal  strip  currents  have  the  same  magnitude  and  only 
differ  in  phase.  This  result  holds  for  both  the  rectangular  and  the 
triangular  array  grids. 

The  relative  TEM  current  distribution  of  the  dominant,  plus- 
minus,  mode  is  shown  in  Figs.  11  and  13.  It  is  seen  that  the  maxi- 

O 

mum  deviation  of  the  current  phase  angle  is  about  12  .  As  stated 
previously,  the  distribution  of  the  plus-plus  mode,  which  is  orthogo¬ 
nal  to  the  plus-minus  mode,  is  determined  by  adding  180°  to  the 
values  shown  in  Figs.  11  and  13. 


The  most  interesting  feature  observed  in  Figs.  (5-8  and  12)  is 
that  the  propagating  TM  mode  may  cut  off  with  scan  angle  before  the 
endfire  grating  lobe  appears  in  real  space.  The  possibility  that  this 
cutoff  causes  the  dipole  array  blindness  is  being  investigated. 


APPENDIX:  TEM  MODE  ORTHOGONALITY 


It  will  be  shown  for  the  general  case  of  multimode  strip  current 
expansions  that  the  capacitance  matrix  in  (37)  is  Hermitian,  and  that 
its  eigenvectors  represent  strip  potentials  which  generate  two  orthog¬ 
onal  TEM  fields. 

The  potential  in  the  unit  cell  of  Fig.  3  is  described  by 


•Kx,y)  =  G(x/y/-^,y,;kt0)n1(y')  +  G(x,y,^  .y';kt0)n2(y')  (Al) 


where  G  is  defined  as  in  (25),  Gg  as  in  (13),  and  2  are  t^ie  str*P 
charge  densities. 

Imposing  the  boundary  conditions  on  the  conductors  yields  the 
integral  equations  for  2 


for  ly 1  <  j 


(A2) 


=  *(-.  y), 

Z  2 

where  2  are  the  constant  potential  values  on  the  strips. 

In  solving  (A2),  the  charge  densities  2  are  approximated  by 
the  finite  series. 


ni(y,)  =  pLi  “ipVy,) 


n2(y')  =  =i  a2p^p(y')« 


(A3) 
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where  i|<p(y')  are  real,  linearly  independent,  basis  functions. 

Substituting  (A3)  into  (A2)  and  projecting  both  sides  of  (A2)  on 
to  the  basis  functions  ^(y')  results  in  a  set  of  2P*2P,  linear,  inhomo¬ 
geneous  algebraic  equations, 

^l'V*  =  ,y^'  V* 

(q  =  1.2 . P)  (A4) 

^2'V  = 

2 

where  <f,g>  is  defined  in  (17). 

Equations  (A4)  can  be  written  in  the  partitioned  matrix  form, 
$  =  M  or.  (A5) 

where 


(A6) 

is  Hermitian  due  to  symmetries 

G(x\y,  x',y’)  =  G*(-x',y',  -x\y) 

(A7) 

and 

G(x,y,  x',y')  =  G*(x',y\  x,y). 

(A8) 

The  elements  of  the  submatrices  A  and  &,  in  (A6),  are 

aqp  =  <*q(y).  G(-  J  .  Y.  ‘  ^  .  Y*)  «l»q(Y')> 

(A9) 

bqp  =  <*q(y).  g(-j  .  .  y')  *p(y')>. 

(A10) 

and  the  vectors  a  and  $  are  defined  by 

T 

?  ”  ®2p*  ®2l * ■ ■ 1 '  ®2P^ 


and 


(All) 


» 


$  =  (^l,^ .  ♦1<l.Vp>.  *2<lt'rl>‘ 

Assuming  that  M  is  non-singular, 


a  =  M  1  <£ 


where 


M-X  = 


4  I  § 

i+  I  I 


(A12) 


(A13) 


(A14) 


is  also  Hermitian.  By  using  the  partitions  of  M  ,  (A13)  can  be  split 


®i  =  4  1  Mo  1 1 


a  l  Mo  4  I 


where 


“i  =  (af . aiP} 


(i  =  1.2) 


i  =  (<i,v . <i.fp>)i 

The  charge  on  the  strips  per  unit  length  in  z  is 


r  -p 

Q.  =  I  a.  <1,4*  >  =  y  of- 
1  p=1  IP  '  P  1  _i 


(i  =  1,2). 


(A15) 


(A16) 


(A17) 


(A18) 


Then  premultiplying  (A13)  with  x  yields 

Q1  =  4  I>  +  S  & 

Q2  =  B+  l)  +  4»2^T  4  *)• 

which  can  be  written  in  the  matrix  form, 


(A19) 


©•<) 


(A20) 


where  the  capacitance  matrix,  C,  is  Hermitian  because  ^  is  real. 
Therefore,  the  eigenvalues  of  the  matrix  are  real  and  the  correspond¬ 
ing  eigenvectors  are  mutually  orthogonal.  In  general,  (A20)  has  two 
linearly  independent  solutions,  which  implies  that  two  linearly  inde¬ 
pendent  TEM  modes  exist  in  the  unit  cell  of  Fig.  3. 
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It  will  now  be  shown  that  when  the  strip  potentials  are  chosen 
equal  to  the  eigenvectors  of  C  their  associated  TEM  field  distributions 
are  orthogonal  over  the  unit  cell. 


Let  K  be  the  two  eigenvalues  of  G,  let  (  .  1  and  (  .  )  be  the 

a'b  -  *a2  *b2 

associated  eigenvectors,  let  <t>a  ^Cx.y)  be  the  corresponding  electro¬ 
static  potentials  in  the  unit  cell,  and  let 


I  =  If  V*  (x,y)  •  V<fc(x,y)  dxdy . 

unit  cell  a  D 


(A21) 


It  is  desired  to  show  that  1  =  0.  The  surface  integral  of  (A21) 
can  be  converted,  via  Green's  Theorem,  into  the  line  integral 

a*!!(x,y) 

I  =  $  <J>a(x,y)  — gjj - dx,  (A22) 

G 

where  the  integration  path  C  runs  counterclockwise  over  the  unit  cell 
boundary  and  clockwise  over  a  closed  contour  hugging  the  strips. 
By  Floquet's  Theorem,  the  contributions  to  the  integral  from  opposite 
sides  of  the  unit  cell  boundary  cancel.  Therefore,  (A22)  equals  the 
contribution  from  the  strips,  that  is 

h 
* 

i  =  /  vi  y) 

n 

■  i 

h 

2 

+  J  <P 
h 

'  2 


/W 

i(7 


y) 


H(?"y) 

3X 


a<(^+,y) 

3x 


dy 


(A23) 


a*b("2'  y) 


3X 


3X 


dy 
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(the  normal  points  toward  the  strips).  The  terms  in  the  brackets 
represent  a  discontinuity  in  the  normal  electric  field  across  the 
strips.  This  discontinuity  is  proportional  to  the  strip  charge  den¬ 
sity.  Therefore, 


1  *  r  V>  *  r  <V?  %2>' 


(A24) 


where  an<*  '1^2  are  tbe  resPective,  strip  charge  densities  asso¬ 


ciated  with  eigenvector 


V  ). 


^2 


By  expanding  according  to  (A3), 


one  finds  that 

1  =  t  pi °*blP  <4>a(‘  r  'y)-V  +*7  pi  °b2p  -y)'V'  (A25> 

and,  therefore,  from  (A4) 

♦a-|  P  *  ♦aO  P  * 

I  =  —  I  cr>  ,  <1,(11  >  +  —  l  a,  „  >. 

eo  p=l  blp  eo  p=l  b2P  p 

Substituting  from  (A18), 

*  * 

t  =  *alQbl  .  (|>a2Qb2 

1  ~  Z  £ 

O  O 


(A26) 


(A27) 


which,  from  our  previous  definition  of  the  eigenvalues  and  eigen¬ 
vectors  of  Q,  can  be  written 


*  * 

I  -  +al*b*bl  *a2*b*b2 

e„  cn 

o  o 

However, 

*al*bl  +  *a2*b2  =  0 


(A28) 


since  the  eigenvectors  of  Q  are  orthogonal  and  consequently  1  =  0. 
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PHASE  SHIFT  WALLS 


Fig.  2  Strip-line  model  of  dipole  array  element. 


I - 1 
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Fig.  3  Rectangular -grid  unit-cell  waveguide  with  balanced  strip 

transmission-line  loading. 


Modified  rectangular  unit-cell  for  the  analysis  of  feed  region 
inodes  in  triangular-grid  arrays. 
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SCAN  ANGLE  0  DEGREES 


Phase  constant  of  the  propagating  TMZ  mode  vs.  scan  angle 
for  the  loaded  unit-cell  waveguide  of  Fig.  3,  showing  its 
dependence  on  dy  (H-plane  spacing). 
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SCAN  ANGLE  9  DEGREES 


Fig.  7  Phase  con stant  of  the  propagating  TMZ  mode  vs.  scan  angle 
for  the  loaded  unit-cell  waveguide  of  Fig.  3,  showing  its 
dependence  on  w  (strip  separation). 
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SCAN  ANGLE  9  DEGREES 


Phase  constant  of  the  propagating  TMZ  mode  vs.  scan  angle 
for  the  loaded  unit-cell  waveguide  of  Fig.  3,  showing  its 
dependence  on  h  (strip  width). 
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SCAN  ANGLE  9  DEGREES 


Fig.  9  Phase  constant  of  the  propagating  TMZ  mode  vs.  scan  angle 
for  the  loaded  unit-cell  waveguide  of  Fig.  3,  showing  its 
dependence  on  the  strip  current  distribution. 
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SCAN  ANGLE  Q  DEGREES 


Fig.  11  Phase  difference  of  the  plus -minus  TEM  strip 

currents  vs.  scan  angle  for  the  loaded  unit-cell  waveguide 
of  Fig.  3. 
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Phase  constant  of  the  propagating  TMZ  mode  vs.  scan  angle 
for  the  loaded  unit-cell  waveguide  of  Fig.  4.  (triangular  grid) 
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Fig.  13  Phase  difference  of  the  plus -minus  TEM  strip 

currents  vs.  scan  angle  for  the  ioaded  unit-cell  waveguide 
of  Fig,  4. 


